ОПТИМИЗАЦИЯ СТРОИТЕЛЬНЫХ КОНСТРУКЦИЙ НА ОСНОВЕ ГЕНЕТИЧЕСКОГО АЛГОРИТМА 

А.Г. Юрьев, С. В. Клюев, А.В. Клюев 

Белгородский государственный технологический университет им. В. Г. Шухова 
Е-таіІ: КІуиуеѵ@уапсІех.ги 

Предложена методика оптимального проектирования несущих конструкций на основе генетического алгоритма. В качестве при- 
мера рассмотрено проектирование стальной рамы при варьировании 9 параметров с использованием метода конечных элемен- 
тов. Выявлен наилучший вариант, соответствующий минимуму объема материала рамы. 


В последние десятилетия в области инженерно- 
го дела, а также экономики и планирования наме- 
чается стремительный переход от допустимых ин- 
женерных и управленческих решений к оптималь- 
ным решениям. Однако современная теория опти- 
мизации пока не удовлетворяет требованиям инже- 
нера-проектировщика в связи с тем, что ее строгие 
математические методы не учитывают реальных 
ситуаций проектно-конструкторских задач. Вместе 
с тем современная, все более усложняющаяся 
практика проектирования и конструирования нуж- 
дается в эффективных математических средствах 
решения таких задач. 

Характерной чертой нового подхода является 
комплексная разработка, позволяющая проектиро- 
вать систему в целом, а не по отдельным ее частям. 
Поэтому одной из чрезвычайно важных научных и 


прикладных задач является разработка методоло- 
гии оптимального проектирования сложных тех- 
нических систем - системного проектирования. 

Конструкцию характеризует ряд показателей: 
стоимость, надежность, вес, габариты, время раз- 
работки и др., которые могут находиться во взаим- 
ном противоречии. Трудность решения задачи со- 
стоит в недостатке априорной информации, 
необходимой для поиска оптимального варианта 
конструкции. Поэтому процедуру проектирования 
целесообразно строить так, чтобы на каждом его 
последующем этапе объем информации о кон- 
струкции возрастал. В то же время необходимо ис- 
ключать неудовлетворительные варианты, выяв- 
ленные в ходе проектирования. Таким образом, 
должны сочетаться две тенденции - генерация 
многообразия вариантов и усечение выявленного 
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множества [1]. Предлагаемая процедура проекти- 
рования согласуется с эволюционной стратегией 
оптимизации и, в частности, генетическими алго- 
ритмами (ГА) [3]. 

Проектирование конструкции представляется в 
виде некоторой последовательности уровней ее 
разработки, которые характеризуются степенью 
детализации ее элементов. Такой способ проекти- 
рования можно связать с некоторой иерархической 
моделью, обладающей спецификой иного класса 
структурных эволюционных моделей - деревьев 
последовательных решений. 

Генетические алгоритмы получили развитие в 
середине 60-х гг. XX в. благодаря трудам Дж. Хол- 
ланда. Они имитируют эволюционный процесс при 
особом акценте на генетические механизмы, т. е. 
наследование и рекомбинацию генов. Осуществля- 
ется это некоторым числом (популяцией) искус- 
ственных хромосом (индивидов). Каждая хромосо- 
ма содержит п генов, которые соответствуют п 
искомым переменным оптимизационной задачи. 

Генетические алгоритмы, как и эволюционные 
алгоритмы вообще, применяют для поиска гло- 
бального экстремума функции многих перемен- 
ных. Принцип их работы основан на моделирова- 
нии некоторых механизмов популяционной гене- 
тики, манипулировании набором хромосом при 
формировании генотипа новой биологической 
особи путем наследования участков хромосом ро- 
дителей, случайным изменением генотипа, извест- 
ного в природе как мутация. 
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Рис. 1. Г енетический алгоритм 


На рис. 1 представлен генетический алгоритм. 
Главная идея эволюции состоит в том, чтобы улуч- 
шать фитнес индивидов первой генерации популя- 
ции, пока не будет достигнут критерий прекращения. 

1. Начальная популяция. Вначале необходимо соз- 
дать начальную популяцию индивидов. Так как о 
целевой функции еще ничего не известно, при- 
мем гены индивидов в ранее установленной 
области случайными и равномерно распределен- 
ными, сопоставимо с методом Монте-Карло. 

2. Оценка индивидов. Каждому из новорожденных 
индивидов необходимо на основании целевой 
функции определить фитнес. После этого мож- 
но начать петлю генераций в стремлении улуч- 
шить фитнес индивидов. 

3 . Селекция - первый шаг к этому улучшению , при 
котором индивиды, по стратегии, выбираются 
случайно или на основе их прежнего фитнеса. 
Этот выбор служит либо продолжением рода, 
либо передачей элитного статуса. Индивид с 
элитным статусом нельзя ни исключить из фак- 
тической генерации, ни изменить. 

4. Воспроизводство означает продолжение рода 
отобранных для этого индивидов. Простейший 
способ размножения - копирование индиви- 
дов. Стратегия ГА пытается помимо этого путем 
рекомбинации генов на хромосомах генериро- 
вать лучшие индивиды. Как правило, при этом 
из родительских индивидов образуются пары, 
которые обмениваются между собой генами 
случайным образом и образуют таким образом 
два новорожденных. 

5. Мутация изменяет отдельные гены генопула. 
Благодаря этому при ГА генопул должен обно- 
вляться, потому что с размножением генов шла 
бы потеря многообразия генопула обычно уже 
после малого числа генераций. 

6. Оценка индивидов. После того, как путем реком- 
бинации и мутации хромосомы индивидов из- 
менены, необходимо для каждого из новорож- 
денных определить фитнес. 

7. Замена индивидов. В конце петли генерации 
необходимо установить, какие индивиды из по- 
пуляции исключаются. Иначе без замены попу- 
ляция росла бы все дальше. ГА заменяют обыч- 
но большинство родителей новорожденными. 

8 . Критерий прекращения определяет длительность 
процесса оптимизации и имеет решающее зна- 
чение для оценки результата. Используются два 
варианта формулировки критерия прекраще- 
ния: 1) контроль фитнеса, в результате которого 
процесс прекращается, если в рамках устано- 
вленного числа генераций максимальная вели- 
чина фитнеса в популяции кардинально не 
улучшается; 2) установление числа генераций. 
Оба варианта имеют несовершенства. При пер- 
вом варианте может случиться, что максимальная 
величина фитнеса в популяции долгое время не из- 


62 








Математика и механика. Физика 


меняется, а затем, например, вследствие удачной 
мутации, наступает ее улучшение. Поэтому в слу- 
чае слишком раннего прекращения достигается 
лишь субоптимум. 

При втором варианте критерия оптимум вооб- 
ще не ставится во главу угла. Часто выбирается ре- 
шение, удовлетворяющее в достаточной мере допу- 
скам, установленным для варьируемых параметров. 
В этом случае все же необходимо проведение опти- 
мизационных процессов с различным числом гене- 
раций, чтобы можно было дать определенную 
оценку результата. 

Путем селекции и исключения индивидов с от- 
носительно плохим фитнесом популяция также пе- 
реоценивается. Разумеется, вследствие этого гено- 
пул может потерять отдельные хорошие гены. ГА 
пытаются сделать эти потери по возможности ма- 
лыми [1,2]. 

В качестве примера рассмотрим проектирова- 
ние стальной рамы. Минимизируется объем, а вме- 
сте с этим вес и расход материала. Задаются допол- 
нительные условия, касающиеся условий прочно- 
сти и устойчивости, общей и локальной. В качестве 
материала применяется сталь С345 со следующими 
характеристиками: плотность р=78,5 кН/м 3 ; мо- 
дуль продольной упругости Е= 2, Г ІО 5 МПа, коэф- 
фициент поперечной деформации ѵ=0,3; предел 
текучести К у = 360 МПа, расчетное сопротивление 
д-300 МПа. 

Заданы размеры рамы: (=4,8 м, 1=4 м и интен- 
сивность нагрузки на горизонтальную проекцию 
ригеля 8 кН/м. Оптимизируется 9 параметров си- 
стемы (рис. 2): Ь - половина ширины полки дву- 
тавра; к щ - высота сечения ригеля слева; /г рп - вы- 
сота сечения ригеля справа; /г св - высота сечения 
стойки вверху; к т - высота сечения стойки внизу; 
? рп и / рс , ( п и ( с - толщина полки и стенки ригеля и 
стойки соответственно. 
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Рис. 2. Профиль рамы и сечения ригеля и стойки 


Из условий соблюдения локальной устойчивости 
элементов двутаврового сечения установлены допу- 
скаемые отношения их размеров: для полки - 10, для 
стойки - ПО. Верхняя и нижняя границы искомых 


переменных устанавливаются следующими (в см): 
5<6<15; 10</г рл <30; 10</г рп <30; 0,5<(„<1,5; 0,5<( с <1,5; 
10</г св <30; 10</г сн <30; 0,5</ сп <1,5; 0,5</ с <1,5. 

Целевая функция, выражающая объем материа- 
ла рамы, имеет вид 

^ = 2, 4 ? сс (й св + /? сн - 2 Г сп ) + 

+2 (% л + V “ 2 ? рп ) + Ь(\ 9, 2 ( сп + 16 Г рп ). 

Это же выражение с обратным знаком опреде- 
ляет фитнес-функцию. 

Используя для расчета метод конечных элемен- 
тов, разбиваем половину объема рамы на 108 элемен- 
тов - 48 в ригеле и 60 в стойке, при этом в попереч- 
ном сечении - 6 элементов (по 2 в полке и стенке). 


В соответствии с числом оптимизируемых пара- 
метров вводим модель хромосомы 
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Элементы сведены в 4 группы - полка и стенка 
ригеля, полка и стенка стойки. Узлы сведены в 26 
групп, передвижения которых соответствуют генам 
хромосомы. 


Таблица 1. Размеры элементов профилей 



ъ 

К 

к 

4п 

4с 

К 

К 

4п 

4с 

1 

Б5 

138 

143 

6,4 

5,0 

100 

105 

5,0 

5,0 

2 

69 

209 

231 

5,7 

5,1 

298 

281 

6,6 

5,0 

3 

50 

120 

157 

5,6 

5,0 

102 

158 

5,0 

5,0 

4 

50 

120 

194 

5,6 

5,0 

101 

108 

5,0 

5,0 

Б 

63 

191 

252 

5,2 

5,0 

191 

300 

5,8 

5,0 

6 

50 

115 

123 

7,0 

5,0 

102 

119 

5,0 

5,0 

7 

50 

161 

162 

5,5 

5,0 

101 

121 

5,0 

5,0 

8 

51 

128 

169 

5,6 

5,0 

103 

116 

5,0 

5,1 

9 

51 

131 

154 

6,1 

5,1 

102 

120 

5,0 

5,0 

10 

50 

126 

161 

5,8 

5,0 

110 

117 

5,1 

5,1 


Для расчета рамы применялись следующие пара- 
метры: индивидов 20; генераций 300; селекционное 
давление 0,1; норма мутации 0,05; число замен 15. 

Отклонения для мутаций носили линейный ха- 
рактер и составляли для толщин полок и стенок 
0, 001. ..0, 0001, а для других элементов 0,01. ..0,001. 

Достигнутые результаты (размеры в мм) пред- 
ставлены в табл. 1. В табл. 2 представлены объем 
рамы V, м 3 , наибольшее нормальное напряжение сг, 
МПа, и отношение Ь/( рп и Ъ/і т . 


Таблица 2. Геометрические характеристики рамы и напряжения 



V 

сг 

Ь/1р п 

Ь/1 т 

1 

0,0150 

161 

7,4 

9,4 

2 

0,0265 

160 

11,7 

10,1 

3 

0,0150 

ібі 

8,5 

9,5 

4 

0,0147 

160 

8,4 

9,5 

5 

0,0224 

388 

11,5 

10,4 

6 

0,0153 

338 

6,8 

9,5 

7 

0,0149 

162 

8,6 

9,5 

8 

0,0149 

162 

8,7 

9,7 

9 

0,0152 

162 

7,9 

9,7 

10 

0,0149 

161 

8,2 

9,3 
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Из табл. 2 видно, что отношения Ъ/І щ и Ь/і с „ пре- 
вышены во 2-й и 5-й строках. В 5-й и 6-й строках 
превышено расчетное сопротивление. Таким обра- 
зом, тремя вариантами решения следует прене- 
бречь. 

Толщины стенок соответствуют принятой ни- 
жней границе. Остальные размеры находятся в 
пределах установленных границ. В семи оставлен- 
ных вариантах решения объем отличается незначи- 
тельно (в пределах 3,4 %). Минимальному объему 
соответствует 4-й вариант решения. 

Выводы 

1 . Генетические алгоритмы являются мощным по- 
исковым средством. Решение, полученное на их 
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основе, субоптимально, но это не мешает при- 
менять алгоритмы для поиска глобальных эк- 
стремумов при оптимизации строительных 
конструкций. 

2. Генетические алгоритмы по сравнению с из- 
вестными аналитическими методами оптимиза- 
ции приемлемы для решения многопараметри- 
ческих невыпуклых задач. 

3. Решение можно уточнить, увеличив густоту сет- 
ки конечных элементов. При увеличении коли- 
чества оптимизируемых параметров возрастает 
число индивидов и генераций. Вместе с тем по- 
вышаются затраты машинного времени, что в 
ряде случаев может служить показателем оцен- 
ки использования генетических алгоритмов. 
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